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Abstract 

Numerous observational studies have revealed the ubiquitous presence of multiple stellar populations in glob¬ 
ular clusters and cast many hard challenges for the study of the formation and dynamical history of these stellar 
systems. In this Letter we present the results of a study of the kinematic properties of multiple populations in 
NGC 2808 based on high-precision Hubble Space Telescope proper-motion measurements. In a recent study, 
Milone et al. have identified five distinct populations (A, B, C, D, and E) in NGC 2808. Populations D and E 
coincide with the helium-enhanced populations in the middle and the blue main sequences (mMS and bMS) 
previously discovered by Piotto et al.; populations A, B, and C correspond to the redder main sequence (rMS) 
that in the Piotto et al. was associated with the primordial stellar population. Our analysis shows that, in the 
outermost regions probed (between about 1.5 and 2 times the cluster half-light radius), the velocity distribution 
of populations D and E is radially anisotropic (the deviation from an isotropic distribution is significant at the 
~3.5-<r level). Stars of populations D and E have a smaller tangential velocity dispersion than those of pop¬ 
ulations A, B, and C, while no significant differences are found in the radial-velocity dispersion. We present 
the results of a numerical simulation showing that the observed differences between the kinematics of these 
stellar populations are consistent with the expected kinematic fingerprint of the diffusion towards the cluster 
outer regions of stellar populations initially more centrally concentrated. 

Subject headings: proper motions — stars: population II — (Galaxy:) globular clusters: individual (NGC 2808) 

— Galaxy: kinematics and dynamics 


1. INTRODUCTION 

Photometric and spectroscopic studies over the last 15 years 
have revolutionized our understanding of globular clusters 
(GCs). Once thought to be “simple stellar populations” with 
a single age and composition, essentially all GCs host mul¬ 
tiple stellar populations (MSPs), as revealed by their multi¬ 
ple photometric sequences. The exquisite precision and sta¬ 
bility of the Hubble Space Telescope (HST ) has made much 
of th is possible (s ee Milone et all 12015 a b: Na rdiello et all 
120151 iPiotto et all 120151 and references therein). MSPs 
can be traced to different compositions in terms of light 
elements (such as Na, O, Al, Mg observed spectroscopi- 
cal lv: see iGratton et all 120121 and references therein) and 
He iPasauini et al.l 1201 ll IVillanova et al.l 120121 iMarino et all 
120141) . These observational findings have cast a number of 
challenges for the study of the formation and evolution of 
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GCs. 

In order to make progress toward a complete picture of the 
formation and evolutionary history of GCs, as well as to con¬ 
strain the possible paths for the theoretical study of these com¬ 
plex challenges and questions, a major effort to combine ob¬ 
servational data is essential. Spectroscopic and photometric 
studies have begun to shed light on the chemical properties, 
number of distinct sequences and range of ages of different 
stellar populations. 

Structural and kinematic properties of different stellar pop¬ 
ulations are two additional key pieces of the puzzle, as they 
contain essential information to build a complete picture of 
MSP cluster formation and dynamical evolution. Accord¬ 
ing to a number of d i fferent formation scen a rios (see, e.g„ 
D Ercole et al.l 12008 1: iDecressin et all 120071 : iBastian et al. I 
20131) . second-generation (2G) populations should form more 
conce ntrated in the cluster inner region^] Numerical simula¬ 
tions (IVesperini et al.l[2013l) have shown that several clusters 
should still retain some memory of this initial spatial seg¬ 
regation and a few observational studies have indeed found 
clusters in which 2G stars are still concentrated in the clus¬ 
ter inner regions dSqlh martakl [2007t |BeHjni_etalJ |2009a|; 


Johnson & Pilachowskil2012F Milone et all2 012tlLardo et all 

201 lHBellini et alJl2013aHCordero et al 1120 1 41) . 


Little is however known to-date on the proper-motion (PM) 
based kinem atics o f MSPs (wi th the only exc eptio ns of u> Cen 
and 47 Tuc: iBellini et al Jl2009bl lAnderson & van der Marell 
120101: iRicher et al.l 120131) . Richer et al. (2013) probed the 
outer regions of 47 Tuc (at ~ 1.9 it) and found that 2G stars 
are characterized by a radially anisotropic velocity distribu¬ 
tion while the 1G population is isotropic. As discussed in Sec¬ 
tion 0] this trend is consistent with the expected kinematical 


1 Hereafter, we will use “populations” and “generations” as synonyms. 












































2 


implications of the formation models mentioned above; the 
study presented in this Letter is aimed at exploring whether 
the kinematics of MSPs in NGC 2808 further confirms this 
general expectation and the trend identified in 47 Tuc. 

NGC 2808 is one of the most massive Galactic GCs and one 
of the most massive Galactic GCs and one of the few known to 
host a super-He-rich subp opulation (see, e.g jD’Antona et all 
120051 IPiotto et ahl 120071) . The initial pictu re of three dis¬ 
tinct stellar populations found bv lPiotto et al.l (120071) has been 
recently shown to be even more complex bv iMilone et al.l 
( 2015b ). who identified five distinct populations (named A, 
B, C, D and E) along the red-giant branch (RGB) and the 
main sequence (MS) of NGC 2808. We point out here that 
populations D and E correspond to the 2G populations iden¬ 
tified by IPiotto et ahl (120071 their mMS and b MS). Popula¬ 
tions A , B, and C correspond to the rMS of the IPiotto et alJ 
( 2007 ) study (that they assumed to be the first generation, 1G). 
New WFC3/UVIS UV data (GO-12605, PI: Piotto) allowed 
us to further split the rMS into the three populations A, B 
and C. Available spectroscopic data discussed in lMilone et alJ 
d2015bl) suggest that population B has Na and O abundances 
usually associated with 1G stars; population C includes stars 
with slightly enhanced Na and slightly depleted O abun¬ 
dances, and is photometrically distinct from population B 
along the RGB tsee IMilone et al.l [2015bi their Fig. 8). No 
spectroscopic data are available to further characterize popu¬ 
lation A. 

In this Letter we shed further light on the properties of this 
complex cluster and present an analysis of the internal PMs 
of its five stellar populations. The outline of this Letter is as 
follows: in Section |2] we describe the data sets and reduction; 
in Section[3]we present our results on the kinematic properties 
of the five populations identified in NGC 2808. In Section [4] 
we discuss the interpretation of the observational findings. We 
summarize our conclusions in Section [5] 

2. DATA SETS AND REDUCTION 

We have compiled high-precision PM catalogs for the cen¬ 
tral regions of 22 GCs, mostly based on the _archival HST 
proposal AR-12845 (P.I.: Bellini, see [Bellini et al.l [2014 ). 
PMs are obtained follow ing the central overlap method 
(lEichhorn & Jeffervslll971 ). via a careful data-rejection pro¬ 
cess. Extensive simulations have demonstrated the reliabil¬ 
ity of the estimated PMs and their errors. The observations 
used to comp ute PMs in NGC 2808 are listed in Table 11 of 
iBellini et al.l d20l4l) . We additionally included WFC3/UVIS 
observations of GO-12605 (6 x 650s in F336W and 6 x 97 s 
in F438W), taken in 2013, to further extend the available time 
baseline and increase the PM accuracy. The catalog con¬ 
tains over 86 000 stars down to ~ 5 mag below the turn off 
(TO). Well-measured stars have a typical PM error of less than 
~ 3 0 /ias yr" 1 (~1.3 kms -1 , at a distance of 9.6 kpc, lHarrisI 
1996). 

The photometric catalog of IMilone et al.l ( 12015b ) com¬ 
bines the GO-10775 F6 06W and F814W ACS /WFC pho¬ 
tometry presented in I Anderson et al.l ( 120081) with the 
WFC3/UVIS F275W, F336W and F438W exposures of 
GO-126 05. All observa t ions w ere corrected for CTE 
effects (lAnderson & Bedinl 120101) and reduc e d usi ng the 
software family ima2xvm (lAnderson et all 12006) . em¬ 
ploy] n g spat i a 11 y - vary i n g, time-dependent empirical PSFs 
(IBellini et al.ll2013~bl) . Stellar positions were corrected for ge¬ 
ometric di stortion following the recipes in Bellini & Bedinl 
d2009h and lBellini, Anderson & Bedinl (1201II) . Photometry is 


ID VELOCITY DISPERSION 



Figure 1 . (a) ID vs. r for the 5 distinct populations. We least-squares 
fitted a 3rd-order polynomial to the population B profile ((cr®)) (black curve). 
The half_light radius, 1.5 x r^, and 2 x rh (dashed and dotted lines) are also 
indicated. Panels (b) to (e) show the normalized difference between a^ of 
populations A, C, D, E, respectively, with respect to that of the reference 
population B. Colored points refer to bins containing the same number of 
stars per population. Black points refer to 3 fixed radial intervals: (1) r < r h, 
(2) rh < r < 1.5 x rh, and (3) 1.5 x rh < r < 2 x rh. See the text for details. 


calibrated as i nlBedin et ahl ( 2005 ). 

We refer to IBellini et al.l ( 2014 ) and IMilone et al.l (1201 5bl) 
for a detailed description of the data reduction. 

The cross-identification of common stars between the PM 
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RADIAL VELOCITY DISPERSION 



TANGENTIAL VELOCITY DISPERSION 
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Figure 2. Similar to Fig.[T]but for the radial component of the motion. See Figure 3. Similar to Fig.[T]but for the tangential component of the motion, 
the text for details. See the text for details. 


and the photometric catalog left us with near 38000 objects 
down to about 3 mag below the TO. In order to analyze in 
detail the internal kinematics of NGC 2808 stars, particular 
care must be taken in removing any PM measurements af¬ 
fected by systematic effects. To this aim, we closely followed 
the recipes described in detail in Section 7.5 of iBellini et al.l 
(120141) . Our final sample consists of about 27000 stars within 
115" (the outermost corner of the field of view) with high- 
precision, high-quality measurements. 


3. THE INTERNAL KINEMATICS OF THE FIVE SUBPOPULATIONS 

To avoid edge effects, we consider hereafter only stars 
within 2 x rh = 96" (Harris 1996,2010 edition). Due to small- 
number statistics, the kinematics of the MSPs on the RGB 
cannot be studied with sufficient detail. On the other hand, 
we have plenty of MS stars at our disposal, in particular: 
598 population-A stars, 3071 B, 1875 C, 4485 D, and 1814 
E. W e will use the same nomenclature and color-coding of 
iMilone et al.l (l2015bh . 
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The ID average transverse velocity-dispersion profile a M of 
each population, obtained by combining together fi a cos S and 
Us, is shown in panel (a) of Fig. [I] We binned each population 
in such a way as to have the same number of measurements (> 
150) in each bin. We marked the position of r h (inner dashed 
line), 1.5 x rh (dotted line), and 2 x (outer dashed line). All 
errorbars refer to l-cr errors. 

The five profiles largely follow the same trend. To better ex¬ 
amine small velocity-dispersion differences among them we 
proceeded as follows. We chose population B (He-normal, 
1G stars) as a reference, and we least-squares fitted a 3rd- 
order polynomial to its ct m profile (red points in panel a) 
to obtain the average trend (<r®). Then, we computed the 
differences Act* between the a t , of the other 4 populations 
(X = A,C,D,E) and the average trend of population B at the 
same radial distance. These quantities are then normalized to 
the average trend of population B itself. These normalized 
differences are shown in panels (b) to (e), as a function of r, 
for populations A, C, D and E, respectively. 

In addition, we binned each population into 3 radial inter¬ 
vals: (1) r < fh, (2) q, < r < 1.5 x r^, and (3) 1.5 x < r < 
2 x rh. The respective values of Aa*/(a^) for these radial 
bins are shown as black points with errorbars. Populations A 
and C (panels b and c) exhibit a slightly increasing Act^/(ct®) 
values the larger the radial distance, while populations D and 
E (panels d and e) seem to have an opposite behavior. These 
differences are marginally significant (at the 2 -ct level). 

We repeated the same analysis shown in Fig. |T| indepen¬ 
dently for the radial and tangential components of motionQ 
Results are shown in Fig.[2]and[3]for the radial (cr ra d) and tan¬ 
gential (<T t an) velocity-dispersion profiles, respectively. 

First, we note the overall trend of CT tan (Fig- EJi) being 
smaller than cr r a! \ (Fig . Eh) for r > />,, in agreement with what 
found by IWatkins et all ( 2015 ) for the RGB of this cluster. 
The radial velocity dispersions of the 5 populations do not 
show significant differences (with only possible hints, at less 
than ~2 -ct level, of a larger cr ra d of population A in the inner¬ 
most region, a larger er ra d for population C in the outermost 
bin), while for r > 1.5 X rh the tangential velocity dispersion 
of populations D and E (panels d and e of Fig. [5} is signifi¬ 
cantly smaller (-0.054 ± 0.019 for D, and -0.086 ± 0.028 for 
E, i.e., at the 2.8- and 3. l-cr level, respectively) than that of the 
reference population B. No significant difference is instead 
found between the kinematics of population B and that of A 
and C, with only a hint (at less than ~2-er level) of a larger 
tangential velocity for populations A and C in the outermost 
regions (panels b and c of Fig [Hi 

To further explore these findings, we computed the radial 
dependence of the deviation from isotropy (cr tan /<j ra d- 1) of 
each population. Results are shown in Fig. [4] The horizontal 
line at 0 corresponds to an isotropic system. As for Fig. H] 
larger black dots refer to data binned in the 3 radial intervals, 
while colored points are used for data binned by keeping the 
same number of stars per bin. 

Populations B and C have an overall similar trend, and 
tend to very slightly deviate from isotropy outside r h. In the 
outermost regions probed (between 1.5 , and 2r h j popula¬ 
tions D and E are radially anisotropic (cr tan /cr la( i- 1 equal to 
-0.093 ± 0.027 for D and -0.140 ± 0.040 for E: a deviation 
from an isotropic distribution at the ~ 3.4-cr and 3.5-cr level, 

2 We decreased the number of bins while keeping the same number of 
measurements within (at least 150) for populations A, D and E. 
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Figure 4. Deviation from tangential-to-radial isotropy (horizontal line) for 
the 5 populations (color-coded as in Fig.[T|. Vertical lines mark the locations 
of rh, 1.5 X rh, and 2 X rh. 


respectively). There is a hint of population A having an oppo¬ 
site behavior with respect to the others: the outer 2 points are 
fairly consistent with an isotropic system, while the innermost 
point indicates a marginally-significant (at the ~2 -ct level) ra¬ 
dial anisotropy for r < r^. As discussed in the Introduction, 
the nature of population A is still unclear and spectroscopic 
data will be necessary to characterize this population in the 
context of the MSPs of NGC 2808. Additional kinematic data 
would also be necessary to shed further light on the kinematic 
properties of this population and clarify the significance of 
possible kinematic anomalies. 


4. THEORETICAL INTERPRETATION 

The two most statistically-significant results presented in 
the previous sections are: i) populations D and E have a ra¬ 
dially anisotropic velocity distribution (Fig. E} in the outer¬ 
most regions, and ii) this difference in anisotropy is due to the 
fact that populations D and E have a smaller tangential ve¬ 
locity dispersion than that of the other populations (Fig. E}. 
No strong differences are instead found in the radial velocity 
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t/t rh (o) 

Figure 5. Time evolution of o’tan/o'rad — 1 measured for particles between 1.5 
and 2 z*h for the 2G (blue line) and the 1G (red line) populations (where r h 
is the projected half-mass radius). The green line in the inset shows the time 
evolution of the mean of the ratio of the 2G to the 1G tangential velocity 
dispersion, while the black line is for the mean of the ratio of the 2G to the 
1G radial velocity dispersion (both measured between 1.5 and 2 r h). Time 
is normalized to the initial half-mass relaxation time (calculated using the 
projected half-mass radius). 

dispersion of all the populations (see Fig.|2|. 

One of the predictions of MSP cluster-formation models, 
in which asymptotic-giant branch (AGB) stars are the source 
of gas out of which 2G stars form, is that the AGB ejecta 
would collect in the cluster central regions and 2G stars would 
initially b e mo re sratially concentrated than the 1G popula¬ 
tion (tD’Ercole et al.l 12008 ). Other formation models, based 
on r apidly-rotating massive stars and/or massive binaries (see 

e.g. iDecressin et all2007t Ide Mink et al.ll2009l:lBastian et al. I 

120 1 3b do not follow the evolution of the processed gas but as¬ 

sume that massive stars are initially concentrated in the clus¬ 
ter inner regions and so is the gas they release. In any case, 
the subsequent long-term dynamical evolution will gradually 
erase the differences in the spatial distribution of 1G and 2G 
stars. 

Here, we present the results of a N-body simulation car¬ 
ried out to explore the evolution of the kinematical proper¬ 
ties of different stellar populations. We started our simulation 
with 50000 equal-mass particles, with 1G and 2G popula¬ 
tions having the same total mass and both following a iKingl 
(119661) model density profile with central dimensionless po¬ 
tential Wo = 7. We assumed that the 2G population is more 
centrally concentrated, with an initial half-mass radius about 
4.5 times smaller than that of the 1G population. As for the 
kinematical properties, both 1G and 2G systems are initially 
isotropic. The cluster is initially tidally limited and is as¬ 
sumed to move on a circular orbit in the external potential 
of the host galaxy (modeled as a point mass). We point out 
that we are not ruling out the possibility that the formation 
process and the subsequent virialization phase might result 
in 1G and 2G subsystems with initially different anisotropy 
profiles, but here we focus solely on the role of relaxation- 
driven long-term dynamical evolution in establishing kine¬ 
matical differences between different stellar populations. This 


simulation has been run with the GPU-accelerated version of 
the NBODY6 code (1 Aarsethll2003t iNitadori & Aarsethll2012l) 
on the Big Red II supercomputer at Indiana University. 

Observational data reveal that the strongest deviation from 
isotropy occurs in the outermost regions for which data are 
available. Figure [5] shows the time evolution of <7 t an/c ra d - 1 
measured for particles between 1.5 and 2 r h (where r h is the 
projected half-mass radius) for the 2G and the 1G populations, 
where cr rac i and cr tan are, respectively, the projected radial and 
tangential velocity dispersions as measured on a plane coin¬ 
ciding with the cluster orbital plane. This is the outermost 
region probed by our observational data where the difference 
in anisotropy between 1G and 2G populations is strongest. 
The results of our simulation show a trend qualitatively con¬ 
sistent with that found in the observational data: as 2G stars 
diffuse from the inner regions (where they formed and are 
initially segregated), they populate the outer cluster regions 
preferentially on radial orbits as shown by the increasing ra¬ 
dial anisotropy of the 2G population. The trend found in our 
simulation a grees wi th what was found in a recent numeri¬ 
cal study bv fHenault-Brunet et al.l (120151) . The different for¬ 
mation and dynamical history of 1G and 2G stars, which are 
currently located in the same outer region of the cluster (in 
this case between 1.5 and 2 rh), is revealed by the differences 
in the kinematic properties. These differences are stronger in 
the cluster’s outer regions, which are initially dominated by 
1G stars. As the cluster evolves, the cluster’s outer regions 
are populated also by 2G stars, diffused from their initial in¬ 
ner location. 

The inset of Fig. 0 shows the time evolution of the mean 
of the ratio of the 1G to the 2G cr ra d and rr tan , measured be¬ 
tween 1.5 and 2rh. The initial differences in the o nll \ and rr uin 
of the two populations evolve in a direction consistent with 
the observed trends: the radial velocity dispersion rr rad of the 
two populations becomes increasingly similar while er tan of 
the 2G becomes smaller than that of the 1G and, in agree¬ 
ment with our observational findings. It is this difference that 
is responsible for the differences in the anisotropy of the two 
populations. 

5. CONCLUSIONS 

In this Letter we have presented an HST study of the kine¬ 
matical properties of the MSPs of_NGC_2808_bascd on high- 
precis ion PM measurements dBellini et al .112014) . In a recent 
study, iMilone et al.l (1201 5 bl) identified five different popula¬ 
tions in this cluster: two of these populations (named D and 
E) correspond to the 2G populations in the mM S and bMS 
identified in an earlier study bv lPiotto et al.l (12007b while the 
other three populations (A, B, and C) correspond to the rMS 
population in the Piotto et al.’s study. In this Letter we show 
that the five stellar populations of this cluster are character¬ 
ized by differences in their kinematics. 

Specifically, we find that in the outermost regions probed 
(between 1.5 and 2 r^) the velocity distribution of the 2G pop¬ 
ulations D and E is radially anisotropic (the deviation from 
isotropy being significant at the ~3.5-cr level). Our data show 
that the larger radial anisotropy in the 2G populations D and 
E is due to the fact that these populations have smaller tan¬ 
gential velocity dispersions than the other populations. Qual¬ 
itatively similar trends have been found bv lRicher et al.ll20l3l 
in 47 Tuc. No significant differences are found in the radial 
velocity dispersions of all the populations and between the 
kinematics of population B and that of A and C (with only 
possible hints of a larger <r ra d of population A in the inner- 
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most region, a larger cr ra( j for population C in the outermost 
bin, and a larger <7 tan for populations A and C in the outermost 
regions all at less than 2-cr level). 

The results of an N-body simulation show that the differ¬ 
ences in both the radial anisotropy and in the tangential veloc¬ 
ity dispersion observed for populations D and E are consistent 
with being the kinematical fingerprint of the diffusion of 2G 
populations from the innermost regions, where they are ini¬ 
tially concentrated, towards the cluster outer regions. A more 
extended survey of models, probing broader ranges of differ¬ 
ent initial conditions, is necessary to fully explore the possible 
evolutionary paths of the 1G and 2G kinematical properties 
and their dependence on initial conditions. 

On the observational side, we emphasize the importance 
of extending the current study to even more external clus¬ 
ter regions. In fact, the outermost regions are essential for 
a complete characterization of a cluster’s kinematical proper¬ 
ties and, in particular, to shed light on the possible effects of 
the external Galactic tidal field. These effects are expected to 
limit the outer development of the radial anisotropy and make 
these more isotropic (and possibly even slightly tangentially 
anisotropic). 

We plan to extend the analysis presented here to other GCs 
as possible, using the PM catalogs of lBellini et al.l (12 0141) in 
combination with the UV photometry of lPiotto et al.l J2015I) . 
This will create a substantial case history that will help us 
in better understanding formation and evolution of MSPs in 
GCs. 

Support for this work comes from STScI grants for HST 
programs AR-12845 and GO-13297. E.V. acknowledges sup¬ 
port also by grant NASA-NNX13AF45G. G.P., S.C., F.D’A. 
and A.R. acknowledge support from PRIN-INAF 2014 (PI: S. 
Cassisi). 
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